Saddle-splay term induced orientational instability in nematic 
liquid crystal cells and director fluctuations at substrates 

A.D. Kiselev^'B 

^Chernigov State Technological University, 
Shevchenko Street 95, 14027 Chernigov, Ukraine 
(Dated: February 2, 2008) 

Abstract 

We analyze stability of the planar orientational structure in a nematic liquid crystal (NLC) 
cell with planar anchoring conditions at both substrates. Specifically, we study the instabilities 
of the ground state caused by surface elasticity at large saddle-splay elastic constant, K24. We 
show that such instabilities are induced by the director fluctuations at confining walls and derive 
the surface part of static correlation functions as a functional integral over these fluctuations 
characterized by the surface free energy. From the surface part of the correlator we derive the 
stability conditions for the planar structure with respect to the fluctuation modes characterized 
by the in-plane wavenumbers and by the parity symmetry. These conditions are analyzed in 
the cell thickness - fluctuation wavelength plane by using the parameterization for the boundary 
curve of the instability region. For relatively small K24, violating the Ericksen inequalities < 
K24 < 2 min(i^i, i^2) the theory predicts that the critical fluctuation mode of the wavelength, Ac, 
will render the structure unstable when the thickness of the cell is below its critical value, dc- 
The parity of the critical mode changes as the twist-splay ratio K2/K1 is passing through the 
unity. Further increase of K24: beyond the second threshold value, 4:KiK2/ {Ki + K2)., leads to the 
instability with respect to the short wavelength fluctuations regardless of the cell thickness. We 
compute the critical thickness and the critical wavelength as functions of -R'24) the twist-splay ratio 
and the azimuthal anchoring strength. 

PACS numbers: 61.30.Gd, 64.70.Md, 61.30.Hn 



I. INTRODUCTION 

Nematic liquid crystals (NLC) confined in 
restricted geometries are technologically im- 
portant Ipj] and have been the subject of in- 
tense studies over the past few decades 0, 01 • 
Anisotropy of the vast majority of NLCs is lo- 
cally uniaxial and molecules of a NLC align 
on average along a local unit director. Ori- 
entational structures in NLCs are thus de- 
fined by distributions of the director n(r) and 
the well-established continuum elastic theory 
provides the phenomenological description of 
orientational distortions V 



In the absence of external fields, orienta- 
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tional structures in spatially bounded NLCs 
crucially depend on the conditions at confin- 
ing walls. These are macroscopically char- 
acterized by the surface contribution to the 
elastic free energy, F,, that adds to the Frank 
elastic energy, F^, describing elasticity of 
NLC in the bulk, to yield the total elastic 
free energy of a NLC in the presence of the 
confining surfaces 



F[n] =F,[n] + F,[n] 



(1) 



{V,nY + K2 (n,V X n)' 



+ Ks[nx (y X n 



df , 



(2) 



1 



iy(n)-i^24[hn)(V,n) 



- {ly, (n, V)n)] 



ds. 



(3) 



where Ki, K2 and are the splay, twist 
and bend elastic constants, respectively; K24 
is the saddle-splay elastic constant; u is the 
outer normal to the surface S; W{n) is the 
surface density of the anchoring energy. 

An important point is that, in addition to 
the anchoring energy which is the anisotropic 
part of the surface tension, there is also the 
elastic contribution to the surface free energy 
that, originally, has been indicated part 
of the elastic energy having the form of a di- 
vergence 0, 0, 1^- This contribution — the 
so-called saddle-splay term (the K24 term) — 
can generally be viewed as the tangential di- 
rector gradient dependent elastic part of the 
surface energy • The other surface elas- 

tic term known as the Ki^ term will not be 
considered in this paper as it can be ignored 
in cases where spatial variations of the den- 
sity and the scalar order parameter are of mi- 
nor importance 11, ^ 3 • 



In the last years K24 specific issues have 
attracted much less attention than the funda- 
mental difficulties caused by the K13 term. In 
particular, though the exact measurements 
of K24 are still missing it was experimen- 
tally estimated to be of the order of the 
Frank elastic constants [l^, EE 113 ■ C)ne 
of the most important theoretical results is 
that the K24 term may induce spontaneous 
twist deformations in hybrid nematic films 
with azimuthally degenerate anchoring con- 
ditions ^3]- Such deformations are man- 
ifested in the formation of periodic stripe 
domains observed in sufficiently thin hybrid 
NLC cells El . 

For planar NLC cells, similar instability 
of the ground state in the presence of the 
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K2A term was considered in Refs. 
Recently, in Ref. 0, Barbero and Perga- 
menshchik suggested that in the proximity 
of the nematic-smectic-A transition the K24 
term grows anomalously large so as to violate 



the Ericksen stability conditions 2^ 1 : 

< K24<2mm{Ki,K2). (4) 

As a result, the uniform equilibrium director 
distribution becomes unstable and a periodi- 
cally modulated nematic phase may occur in 
sufficiently thin planar films. 

The results of Refs. [H H Q are es- 
sentially limited to the special case in which 
the azimuthal anchoring strength is identi- 
cally zero. But the study of possible mecha- 
nisms leading to the formation of modulated 
orientational structures close to the nematic- 
smectic-A transition requires a quantitatively 
accurate description of the instability that 
goes beyond this limitation. At this stage, 
however, even the instability scenario as a 
whole has not been studied in any detail. 

In this paper we intend to fill the gap. Our 
primary goal is the comprehensive study of 
the instability induced by the K24 term in 
the presence of the azimuthal anchoring. 

The idea underlying our general theoreti- 
cal considerations is that instabilities of this 
sort occur when the director fiuctuations at 
confining surfaces become critically diver- 
gent. So, we suggest the method connecting 
the correlation functions of director fiuctua- 
tions and the computational procedure ap- 
plied to perform the stability analysis. This 
method is based on separating out the sur- 
face part of the correlator correlation 
function of the fiuctuation field induced by 
the director fiuctuations at confining walls. 

The layout of the paper is as follows. In 
Sec. ini we express the surface part of the 
static correlation functions of the director 
fiuctuations as a functional integral over fluc- 
tuations at confining walls and explicitly re- 
late the procedure for computing the corre- 
lators to the stability conditions used in our 
stabihty analysis. 

Analytical results for the planar NLC cell 
are described in Sec. IIIII We characterize the 
mirror symmetry properties of the fiuctua- 
tion harmonics and calculate the surface part 
of the correlator. We find that the result is a 
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sum of tfie contributions from the two fluctu- 
ation modes of different symmetry (symmet- 
ric and antisymmetric) and derive the stabil- 
ity conditions for these modes. 

Stability of the uniform planar orienta- 
tional structure is studied in Sec. IIVI We 
analyze the parameterization of the bound- 
ary curve enclosing the instability region 
in the thickness-wavelength plane and show 
that, in addition to the stability interval (jH), 
there are two different intervals for K24: 
(a) 2min(i^i,K2) < K24 < 4K^K2/{K^ + 
K2), where the critical point is character- 
ized by the critical thickness dc and the crit- 
ical fluctuation wavelength Ac; (b) K24 > 
4:KiK2/{Ki + K2), where NIC cells of any 
thickness are unstable with respect to the 
short wavelength fluctuations with A < Aoo- 
It is found that the critical fluctuation mode 
is antisymmetric at K2 < Ki and is sym- 
metric at K2 > Ki. The critical thickness 
and the critical wavelength are computed as 
functions of K24 and the azimuthal anchor- 
ing strength. We also discuss the spectrum 
of director fluctuations at the substrates near 
the critical thickness. 

Finally, in Sec. El we present our results 
and make some concluding remarks. Details 
on some technical results are relegated to Ap- 
pendixes |X] and ^ 

II. CORRELATION FUNCTIONS 
AND STABILITY CRITERIA 

In this section we consider the general pro- 
cedure for computing the correlation func- 
tions (correlators) of NLC director fluctua- 
tions in confined geometry. We remind the 
reader about the standard approach that uses 
functional integrals to represent averaging 
over fluctuations [1, li^. In this approach 
the effect of the confining surface enters the 
theory through the boundary conditions for 
the saddle point equations (Euler-Lagrange 
equations) and for the Green functions. 

We then describe an alternative proce- 
dure, where the part of the fluctuation field 



representing the director fluctuations at the 
surface is separated out by shifting the in- 
tegration variable in the functional integral. 
The corresponding part of the correlator is 
deflned by the surface part of the free en- 
ergy (jni) and involves averaging over the fluc- 
tuations at the surface. Finally, we show that 
the free energy of these fluctuations deter- 
mines stability of orientational structures. 

A. Free energy of fluctuations 

Assuming that the director fleld no deflnes 
the unperturbed orientational structure, we 
shall write the distorted conflguration in the 
form 

n = cos 6 cos no + cos 6 sin (p ni 
-Fsin6'n2, (n^, nj) = 5^, (5) 

where the brackets denote the scalar prod- 
uct. Locally, the vectors no and ni form the 
no-ni plane, so that the angles (p and 6 de- 
scribe in-plane and out-of-plane orientational 
distortions, respectively. 

Substituting Eq. Q into Eqs. © - © will 
give the free energy of the director conflgu- 
ration (0) ClS db functional of the angles and 
9. For small distortions, (j),6 1, Eq. © 
reduces to the familiar form 

n ^ no + 5no, 6no = ni + n2 (6) 

and the lowest order approximation for the 
free energy of director fluctuations — the so- 
called Gaussian approximation — is given by 
the second order term, F^'^\ of the truncated 
series expansion for the energy © 

F[n]^F[no]+F(2)[V.], (7) 
F(2)[V.]=Ff [V^], (8) 

where ip = ^^^^ ~ (^g^ 

component fluctuation fleld; F^"^^ and are 
the bulk and the surface parts of the fluctua- 
tion free energy generated by the correspond- 
ing terms of the free energy ([T]). 
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FIG. 1: Schematic representation of the planar NLC ceh. In-plane and out-of-plane fluctuations 
are described by the angles 4) and 6*, respectively. 



We can now apply the standard varia- 
tional procedure to obtain the saddle point 
equations for the fluctuation free energy. 
(These equations can also be derived as the 
linearized Euler-Lagrange equations for the 
director (0)). In this section, we only need to 
specify the general form of the equations: 



5i)i{Y) 



ir,,Vj(r) = [i^V(r)]. = 0, (9) 



where K is the matrix differential operator. 
Hereafter, matrices and matrix differential 
operators will be indicated by hats and we 
also assume the summation over repeated in- 
dices. 

The general structure of the fluctuation 
free energy can now be expressed as follows 

Jv Js 

(10) 



Js 



where (<^, = ipiAijipj and the surface 
term on the right hand side of Eq. (fTUI) results 
from the integration by parts. 

There are two most important formal 
properties of Fs'^^ and -F/^'*: (a) director 
derivatives along the normal to the surface do 
not enter Fg , and (b) the bulk free energy 



(2) 

Ff^ is represented by the symmetric bilin- 
ear functional Sb, so that = S'^fi/', V'] 
and Sb[^,il)] = Sb[il^,^]. In particular, the 
relation 



[(V',Q(^V)-(¥',Q^'V)]d3, (12) 



which is a version of the Green formula, im- 
mediately follows from the symmetry of Sb- 

B. Correlators and Green functions 

We shall need to write the probability dis- 
tribution of fluctuations at the state of ther- 
mal equilibrium in the form 



P[i/)] = Z-iexp{-/3F(2)[-0]}, 



(13) 



where /3 = [ksT)^^, ks is the Boltzmann 
constant, T is the temperature and Z is the 
partition function given by the functional in- 
tegral 

Z = y"exp{-/3F(2)[^]}I?^, (14) 

where = The components of 

the correlator then be written as the averages 
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FIG. 2: The graphs of the functions 7s(u) (sohd hne) and 'Jaiu) (dashed Hne) in the u-(?24 plane. 
Two cases are illustrated: (a) r < 1 (r = 0.5) and (b) r > 1 (r = 1.5). The endpoint of the 
instability interval for the critical mode u < Umax at q^^ < ^24 < is shown on the left. At 

(2) (3) 

Qc < 924 < Qc , the endpoint of the instability interval for the non-critical mode u > Umin is 
indicated on the right. 



of the fluctuation field products energy is now as follows 



a,(r,rO = [C(r,r')]., = (V^,(r)^,(r')) 

M^)^j{^')P[tP]V'ip. (15) 



= F(2)[V)]-r' / Mr)iJ,i^)dv, (16) 
Jv 

where J = ( "r^ ) ■ relation linking the 

The correlator ()15|) is more appropriately V 2/ 

known as the 2-point static correlation func- correlator and functional derivatives of the 



tion [2^ and so long as we are within the generating functional is 

5^j(r|J 



Gaussian approximation no other correlation C- {r r') = C- (r r'|J)| 

5^ \nZ[3] 



functions are required. 

For functional integrals, explicit analytical 



, , (17) 
j=o 



treatment can be rather involved [27, 28]. For 6J,{r)5Jj{r') j-o 6Jj{r' 

our purposes, it is, however, more convenient where 

to modify the distribution (fT^ by adding a p 

source term to the free energy so as to intro- ^i{r\J) = / '^pi{r)P[^p\3] Vjp. (18) 
duce the generating functional of the correla- 

tion functions through the partition function In our case all the functional integrals are 

of the modified distribution ji^ . The free Gaussian and the saddle point approximation 
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Instability region 




Instability at d < d 



Stability region 



2 3 

r ( = K^/K^ ) 



FIG. 3: Stability diagram in the r~g24 plane. There are two critical values of q24: qc^ (solid line) 
and qc (dashed line). 



yields the exact results. So, the averaged fluc- 
tuation field ()18p can be computed as an ex- 
tremal of the free energy ()16|) which satisfies 
the Euler-Lagrange equations 



i^,,V^,(r|J)=/?-V,(r) 
and meets the boundary conditions 

4iVSi(r|J)L,5 = o, g = gW + Q(^) 



(19) 



(20) 



So long as the boundary value problem 
(Iini)-(|2ni) is linear, the fluctuation field (HH)) 
linearly depend on the source. This relation 
is given by 

^,(r| J) = / ^..(r, r') J,(r') dv\ (21) 
Jv 

where G is the Green function of the operator 
K. Eqs. (|17p and (j2H) show that the correla- 
tor is proportional to the Green function 



C(r,rO=r'G(r,r'). 



(22) 



So, in order compute the correlator C we 
need to solve the boundary value problem for 
the Green function G: 



KikG 



ik^kj [,'■ ) 



r, r 



QikG 



r, r 



\reS 



6{r- r')6ij, 
0, 



(23a) 
(23b) 



where S{r) is the delta function. The prob- 
lem fl23p is at the heart of the conventional 
computational procedures traditionally used, 
e.g., in studies of light scattering in confined 
liquid crystals 2a, M, 21,21,M- 

The key point is that the effects caused 
by the anchoring energy and the surface elas- 
ticity constant are solely incorporated into 
the boundary conditions (j23b|l 



These con- 
ditions affect eigenfunctions (normal fluctu- 
ation modes) and eigenvalues of the opera- 
tor K that can be used to derive the Green 
function in the form of an expansion over the 
eigenfunctions. 

The eigenvalues form the spectrum of fluc- 
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FIG. 4: Stability diagrams in the D-A plane at = Wg and r = K2/K1 = 1.5 for various values 



of 924 = K24/K1: (a) qi^'^ = 2 < = 2.1 < q^/' = 2.4; (b) q]f> = 2.4 < 924 = 2.43 < 2r = 3; and 
(c) 2r = 3 < 924 = 3.3. 



(2) 



tuations and must be positive provided the 
orientational structure ng is stable. Other- 
wise, the functional integrals (fT^ and (fT^ 
do not converge. 



We have thus formulated the spectral 
stability conditions that turn out to be 
closely related to the approach based on 
the generating functional. These conditions 
are in considerable use as stability crite- 
ria in several methods developed to study 
Freedericksz-type transitions in confined liq- 
uid crystals 3J, 35l 



C. 

bility 



Surface part of correlator and sta- 



We now pass on to the approach that 
emphasize the role of the director fluctua- 
tions at the confining surface by using an- 
other transformation of the functional in- 
tegral (jlSp . This transformation has long 
been known as an efficient method to per- 
form Gaussian integrals and we, follow- 
ing the general idea, define the new integra- 
tion variable iff, that vanishes at the surface 
by translating the fiuctuation field i/j: 



(24) 
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2 2.1 2.2 2.3 2.4 

t>24 (K24/K, ) 

FIG. 5: Dimensionless parameter of critical thickness Dc (= dc/le) versus ^24 at r = 1.5 for various 
values of the azimuthal anchoring parameter = W^p/We. 




where 9?^ 



The fluctuation fleld i/? is thus decom- 
posed into the fleld, c^j,, vanishing at the sur- 
face and the fleld that accounts for non- 
vanishing fluctuations at the surface, = 
(fg. When (p additionally satisfles the Euler- 
Lagrange equations Q 



K(p = 0, ip\ =1^1 =ip^ 



(25) 



the fluctuations described by y?^ and ip will 
be statistically independent 



From Eqs. and (fT^ we obtain the 
probability distributions: 

PM = Z,'exp{-f3Fl'\cp,]}, (28) 
PM = Z;'eM-P<!>,[ip,]}, (29) 

where Pb[^b\ is the distribution for the fleld 
of fluctuations in the bulk cp^; and the distri- 
bution characterizes the fluctuations 
at the surface (p^. 

Averaging over the fluctuation flelds (p^ 
and (pfj can now be performed independently. 
As a result, we have 



+ ^] = Ff [(^,] + <|.,[(pj, (26) (AA,) = I A,A,P[il,]Vi, = I A,P,[^,]Vcp, 



(27) 



Derivation of Eq. (j^UI) relies on the formal 
properties listed at the end of Sec. Ill Al It 
involves using the Green formula (|12j) and the 
relation Q^'^VbU ^ 



Vbls=0 



{As)s{A,),, (30) 



where As = As[(Ps] and Af, = Ab[(p^]. 

After substituting the decomposition 
into Eq. |T3|) and using Eq. (^01) to carry out 
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averaging over fluctuations we arrive at the 
expression for tlie correlator in the final form 

ajiry) = Ci^^(r,r') + Clf(ry). (31) 

The two terms on the right hand side of 
Eq. are given by 

Ct^(r,r') = (¥^f(r)v,f(r'))„ (32) 

C^^(r,r') = (¥^.(r|c^j¥.,(rVJ)s, (33) 

where notations for the argument of ip indi- 
cate the dependence of the field (p on cp^ (see 
Eq. (EH)). 

The correlator C^^^ is entirely determined 
by the bulk director fluctuations with the 
probability distribution (j^Hj) . In this case the 
fluctuations at the surface are suppressed and 
the boundary conditions for C^^^ are 

cS\^y)Ls = ^- (34) 



These conditions can be described as the 
limit of infinitely strong anchoring for 
Eq. (l23b|l . The result is that C = C^''^ in 
the case of strong anchoring. 

When the anchoring is not infinitely 
strong, the boundary conditions ()23b|) differ 
from the strong anchoring conditions ()34j). 
Now the fluctuations at the confining wall 
have not been suppressed completely and 
are characterized by the probability distribu- 
tion (123) with the free energy (1771). Ea. (j^ 
shows that these fluctuations transmitted 
into the bulk give rise to the fluctuation field 
(p. Eq. (|HH|) gives the correlator of the field 
(p induced by the fluctuations at the surface 
and determines the difference between C and 
C^^\ In what follows the correlator C^^^ will 
be referred to as the surface part of the cor- 
relator C. 

It is of particular interest that the surface 
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part of the correlator is written as an average 
with the distribution (j29|l . When the free 
energy (j77jl is not positive definite and the 
conditions 



(35) 



are violated, the correlator (jH^ becomes di- 
vergent. So, stability of the orientational 
structure with respect to the director fluctu- 
ations at the confining surface is determined 
by the stability conditions ((23) • Instabilities 
induced by such fluctuations are of our pri- 
mary concern, since neither the anchoring en- 
ergy nor the K2A term may affect the corre- 
lator of the bulk fluctuation field ((221) • 

Technically, using the conditions ((H^ 
enormously simplifies stability analysis as 
compared to exploring the spectrum of fluc- 
tuations js^. The procedure involves two 
steps: (a) solving the boundary value prob- 
lem for the Euler-Lagrange equations ((^: 



and (b) evaluating the free energy of the fluc- 
tuations at the surface ^I7\i . Finally, the sta- 
bility conditions are derived as the conditions 
for the free energy to be positive definite. 

Equivalent computational schemes have 
already been used to study surface elastic- 
ity effects in Refs. HUll We have thus 
developed the method at the very least link- 
ing these schemes and the surface part of the 
correlator that takes into account the direc- 
tor fluctuations at the confining surface. 

In the next section we apply the above 
described procedure to the case of the pla- 
nar orientational structure in NLC cell. As 
an advantage of our approach, the study of 
the fluctuations at the bounding plates will 
be naturally incorporated into the stability 
analysis 
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FIG. 8: Stability diagrams in the D~A plane at = and r = 1.5 for various values of ^24: 



(a) qi^^ = 2 < 924 = 2.K q^/' = 2.4; (b) 
3.3. 



.(2) 



.(2) 



2.4 < ^24 = 2.5 < 2r = 3; and (c) 2r = 3 < ^24 



III. DIRECTOR FLUCTUATIONS IN anchoring potential is 
NEMATIC CELL 



In this section we consider a NLC planar 
cell of the thickness d sandwiched between 
two substrates that are both normal to the z 
axis: z = —d/2 and z = d/2. Anchoring con- 
ditions at both substrates are planar with the 
vector of the easy orientation directed along 
the X axis. Thus, the uniform planar direc- 
tor distribution no = e^, gives the undistorted 
orientational structure. 

The distorted director (0) with rii = e^, 
and n2 = is characterized by the angles (p 
and 6 shown in Fig. ^ For small angles, the 
approximate second order expression for the 



W{n) ^ W{no) + Wee' + W^<P\ (36) 

where Wq and W^i) are the zenithal (po- 
lar) and the azimuthal anchoring strengths. 
Note that Eq. (jHUj) does not imply using the 



where We 



Rapini-Papoular potential |37| . 
W^. It is, however, more reasonable to con- 
sider the case in which the zenithal and the 
azimuthal anchoring energies differ, because 
the surface energy costs for splay-bend and 
twist director deformations in interfacial lay- 
ers are different (recent discussion can be 
found in Ref. Q). 

The cell is invariant under translations 
in the x-y plane and we impose the pe- 
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FIG. 9: Dimensionless parameter of critical thickness Dc (= dc/h) versus at g24 = Qc — 0.1 
for various values of the elatic anisotropy parameter r. Insert at the upper right corner enlarges 
initial decay in the neighborhood of the origin. 



riodic conditions on the angles: il^lx + 
L^,y,z) = il^{x,y,z) and il^{x,y + Ly, z) = 
ilj{x,y,z), where L^. and Ly are the charac- 
teristic lengths of the cell along the x and y 
axes, respectively. We can now write down 
the Fourier series expansion for the fluctua- 
tions: 

^ = 5Z {"^ml^) exp[i{k^x + kyy)] 

m, my>0 

+ '4^*m{z)exp[-i{k^x + kyy)]}, (37) 

where m = {m^jmy), k^ = 2'7rm^/La., ky = 
2TT'my/Ly and an asterisk indicate complex 
conjugation. 

As is shown in Appendix the fluctua- 
tion harmonics with non-vanishing k^ do not 
produce additional instability. So, we shall 
restrict ourselves to the case in which k^ = 0. 
Owing to the translational symmetry, the 
fluctuation harmonics are statistically inde- 
pendent in the Gaussian approximation and 



the free energy of fluctuations takes the form 
of a sum of the energies of the fluctuation 
modes: 

F(2)[V'] = $^(2-5^,o)Fi^)[VJ, (38) 

m>0 

where m = niy. 

Calculation of the free energy is rather 
straightforward (some details and related 
comments are given Appendix ^ and we 
present the result in matrix notations for the 
modified fluctuation harmonics: 

V'™ - = (J °) V'^. (39) 

For brevity, we drop the index m from no- 
tations for the fluctuation field. The energy 
then can be written in the form 

F!^\^] = ^{Sm^]+Sl:^[^]}, (40) 
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FIG. 10: Dimensionless parameter of critical wavenumber Kc (= 2tt/Ac = kch) versus at 

(2) 

(724 = Qc — 0.1 for various values of the elastic anisotropy parameter r. 



u 



r dr [[dr^P]^A + I V^^a.V"^^"'" f = ^'"^ substrates; 

"^""^ 9t- = — — and a cross indicate Hermitian con- 

or 

(41) jugation. The dimensionless parameters used 
in Eqs. dHJ-dlH) are: 



+ -[9>]+fi+'0 + '0^C'V' 



Ai=±l 



A 



r 
1 



C 



T=fJ,U 

1 
r 



5 



1) 



1 
-1 



AiM(g24- (l + r)/2) 



1 

1 



+ 



^w^ 
, We 



(42) 



(44) 



K. 



T 



U 



24 



fc.^, r = -, g24 , 



(46) 
(47) 



(43) 

where Ig and /(i are the zenithal and the az- 
imuthal anchoring extrapolation lengths, re- 
spectively. 



A. Mirror symmetry and parity of 
fluctuations 



In order to evaluate the correlator of di- 
(^45) rector fluctuations and study the stability of 
the planar structure, no = e^;, we shall need 
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FIG. 11: (a)-(b) The critical thickness and (c)-(d) the critical wavelength parameters as a function 
of the twist-splay ratio r for various values of r^. The graphs are computed by using the Cauchy 
relation: 524 = (1 + r)/2. The saddle-splay parameter ^24 is in the instability region (|88p when r 
lies in the intervals: 3 - 2\/2 <r<l/3or3<r<3 + 2^2. 



to solve the Euler-Lagrange equations for the 
fluctuation free energy functional pUjl : 

Lip = 0, (48) 
L = Adl- Bdr - C. (49) 

These equations are invariant under the mir- 
ror symmetry transformation [i21l] : 



^(r)^Pc^(_r), P = 




Algebraically, this result follows because the 
matrices P, A and C are commuting, PA — 
AP = PC — CP = 0, whereas the matrices 
P and B anticommute, PB -\- BP = 0. 

The identity P^ = I, where / is the iden- 
tity matrix, shows that invariant sets of the 



solutions of Eq. (j48|) can be characterized by 
the parity with respect to the transforma- 
tion (jSni): 

V,,,(-r) = ±PV.,aM, (51) 

where ^ind i/j^ will be referred to as 
the symmetric and the antisymmetric fluc- 
tuation modes, respectively. For symmetric 
[antisymmetric] fluctuation fields, the parity 
relation (j51|) means that the in-plane and 
the out-of-plane components, (p ^.nd 6, are 
represented by even [odd] and odd [even] 
functions of r, correspondingly: 0(— t) = 
0(r) [0(-r) = -0(r)] and 9{-t) = -9{t) 

[e{-r) = eir)]. 

The general solution is now a sum of the 
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FIG. 12: Elements of the covariance matrix ((TSJ as functions of u (= kyd/2) for (a) symmetric 
and (b) antisymmetric fluctuations at the lower substrate, z = —d/2. The thickness parameter 
D = Dc + \.^ = 1.27 is noticeably above its critical value Dc = 0.27 at ^24 = 9c^'* + 0.3 = 1.7 < q^f^, 
r = 1.5 and = 1.0. 



symmetric and the antisymmetric modes: 

(^(r) = V'.(r)+Va(r). (52) 

We shall write expressions for the modes in 
matrix notations through the two fundamen- 
tal matrices ^'^'*^(r) and '^^"'\t) composed 
from solutions of the corresponding symme- 
try. These 2x2 matrices satisfy both the 
Euler-Lagrange equation (jiH)) and the par- 



ity relation (j^T|) . In addition, ^^"^(r) will 
be conveniently normalized by the condition: 
^(")(m) = /. It can be verified that solving 
Eq. pH|) yields the following result: 

V^Jr) = ^(°)(r)v?,, a = s,a, (53) 
¥^\t) = ^^iT)[^Mr\ (54) 



Mr) 



cosh T + pr sinh r —pr sinh r 

\ pr cosh r sinh r — pr cosh r 

sinh T + pr cosh r —pr cosh r 

pr sinh r cosh r — pr sinh r 



(55a) 
(55b) 



where y?^ are the vectors of integration con- slants and 

p = (l-r)/(l + r). (56) 
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FIG. 13: Elements of the covariance matrix H75|) as functions of u (= kyd/2) for (a) symmetric 
and (b) antisymmetric fluctuations at the lower substrate, z = —d/2. The thickness parameter 
D = Z?c + 0.02 = 0.29 is in the immediate vicinity of its critical value = 0.27 at 924 = gi^'* + 0.3 = 
1.7 < q'c \ r = 1.5 and = 1.0. 



After substituting Eqs. into 
Eq. ()40j) and using the symmetry rela- 
tion (j3T|) . we find that the symmetric 
and the antisymmetric fluctuations indepen- 
dently contribute to the free energy of the 
fluctuation field (jK^ : 

(57) 

M^ = u{Adr¥^\^^-\B)+Qf. (58) 



We can now combine Eqs. ()43 p - (P3j) and 
Eqs. (IMjl - ffSUj) to derive expressions for the 
matrices Ms and Ma from Eq. (j3Hj) . The re- 
sult is given by 



M. = — 



1 

Js 

1 

Pa 



1 — p)u tanh^ u + W(f,(3s 
u{q24Ps - (1 - p) tanhw) 

(1 - p)m + w^Pa 
u{q2iPa - (1 - p) tanhw) 



u{q24f3s - (1 - p) tanhw) 
(1 - p)u + wePs 

u{q2iPa - (1 - p)tanhn; 
(1 — p)u tanh^ u + Wg(3a 



(59a) 
(59b) 



16 



where 

/5s,a = tanhw =1= p'u(l — tanh^ m). (60) 

B. Director fluctuations at substrates 

Now we apply the analytical results of 
Sec. nil Bl to evaluate the energy of fluctua- 
tion harmonics at the substrates and the sur- 
face part of the correlation function. Accord- 
ing to the procedure described in Sec. Ill CI 
the first step is to solve the boundary value 
problem (j^ . 

To this end we shall find the fluctuation 
field that satisfies the Euler-Lagrange equa- 
tions (jlSj) and the boundary conditions 



(±) 



(61) 



where c^^^-* ['^'•^•*] is the value of the fluctu- 
ation field at the upper [lower] substrate of 
the cell. By analogy with the representation 
given in Eqs. and ()53|). we can write the 
field in the form 

cp{r) = ^W(r)cpW + ¥-\t)cp^-\ (62) 
¥''\uu) = S^J, /i,z/ = ±l, (63) 

where the matrices '^^'^^ and are the so- 
lutions of the Euler-Lagrange equation (|iS|) 
that are equal to the identity matrix at one 
substrate and vanish at the other. 

It can be readily verified that and 
are related to the matrices derived in 
Eqs. fl54|) - ()55b|) through the linear combina- 
tions 

^(+) = + (64a) 

The fluctuation field (jU^ thus has been de- 
fined. 

Clearly, Eq. and Eq. present the 
fluctuation field (p in two different forms. 
This in combination with Eqs. (j64a|) - ()64b|) 



yields the relation between the vectors c^*-^^ 
and (y2 : 



^ - I ^^^^ 



T 



(65) 



where T 



I I 
-P. 

The general expression for the energy of 
fluctuations at the surfaces has been derived 
in Eq. (j2ZI)- In our case this energy is a 
quadratic form of the boundary values of the 
fluctuation field (jU^ . <p^~^^ and (p^~\ that can 
be evaluated by substituting Eq. (jU^ into 
Eq. (j4(Jj) . From the other hand, the energy 
expressed in terms of cp^ and ip^ has already 
been obtained in Eq. (jHZj). So, with the rela- 
tion (jU^ we are left with the following result 



KiS 



T+NT 



Ms 
M„ 



(66) 
(67) 



where the matrix is described in Ap- 
pendix IHl (see Eqs. dnsi-dnsi))- 

After substituting the fluctuation field 
into the expression for the surface part of the 
correlator (jH^ we have 



r, r 



(^i(r)^j-(r')), 



im \ I ran nj \ I 



(68) 



where Dmn = {(p^m(pn^)s- In our case, 
we can easily carry out the averaging over 
fluctuations by calculating Gaussian integrals 
with the probability distribution defined in 
Eqs. (123) and (jnni)- The result is 



D = 



/){++) /){+-) 
£){-+) £){—) 



ksTd 



N-\ (69) 



In Appendix^ the Green function method is 
found to yield the same result. 

By using Eqs. flM|l and (jUTjl the surface 
part of the correlator can now be derived 
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from Eqs. 



in the final form: 



C'^'\t, t') 



ksTd 



2KiS 

f{t,t') = ¥'\t)m;^[¥'\t')] 

+ ¥''\T)M-^[¥''\T')f. 



(70) 



(71) 



As is seen from Eq. ()71|). the surface part 
of the correlator is a sum of two terms that 
represent the contributions coming from the 
symmetric and the antisymmetric fluctua- 
tion harmonics at the substrates. The cor- 
responding terms on the right hand side of 
Eq. (j57p provide expressions for the free en- 
ergies of these harmonics. 

We can further emphasize the role of the 
matrices Mg and Ma- For this purpose, we 
consider the case in which r = r' = ±m. 
Since C^^\±u,±u) = 0, the covariance ma- 
trix of the director fluctuations at the up- 
per and lower substrates is determined by the 
correlator (fTOIl : 



C^'\±u, ±u) = C{±u,±u) 

^ z=±d/2 



(72) 



From Eq. ()7ip we have 



F{u,u) = M:^ + M~\ (73) 
F{-u, -u) = P{M:' + M-i)P. (74) 

The fluctuations at the substrates are thus 
entirely described by the matrices Ms and 
Ma- In Sec. IIV Dl we shall discuss symmetric 
and antisymmetric fluctuations at the lower 
substrate of the cell and present some numer- 
ical results for the elements of the matrix 



(75) 



From Eqs. (|7U|) and (|74j) this matrix is pro- 
portional to the contribution of the corre- 
sponding fluctuation mode to the covariance 
matrix ()72|). 



IV. STABILITY OF PLANAR STRUC- 
TURE 

It is now our task to study the stability 
of the planar structure. In the preceding sec- 
tion we have derived all the analytical results 
required to perform analysis of the stability 
conditions ()35|) efficiently. In our case violat- 
ing these conditions will render the surface 
part of the correlator (7^'') divergent because 
the free energy is not positive definite 
and the corresponding Gaussian integrals di- 
verge. 

Equivalently, Eqs. and (jFFjl show that 
the planar structure is stable only if the ma- 
trices Ms and Ma are both positive definite. 
This yields the following two stability condi- 
tions: 

det M, > and det Ma > 0, (76) 

that determine stability with respect the 
symmetric and the antisymmetric fiuctuation 
modes. 

In order to proceed with our analysis we 
shall need to write down the expressions for 
the determinants. From Eqs. (|59a|) and ()59b|) 

we have 

det Ms = wgWtf) H ^ ^ ^ [wg tanh u + Wfj,) 



-u q2A( q2A - 2 



(3s{u 
(1 — p) tanhu 



(77) 



det Ma = weWfj) + 



u q24 524-2 



(1 — p)u , ,2 \ 

(w^tanh u^Wq) 
1 — p) tanh-u 



(78) 



The determinants are written in the form of 
a sum of three terms, so that only the last 
term can be negative (the functions Ps,a{u) 
are defined in Eq. (p]j) and cannot be nega- 
tive). 

This term is always negative at ^24 < 
and we concentrate on the case of our pri- 
mary concern in which g24 is positive. In this 
case the term will be negative if the inequal- 
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ity 



g24 > 1a{u) = 2 



[1 — p) tanhw 



(79) 



is fulfilled. Given the value of ^24, the values 
of the parameter u that satisfy the instability 
condition (f?^ form the instability interval for 
the corresponding fluctuation harmonics. 

A. General method 

The key point underlying our analysis is 
that the determinant det will be negative 
provided the cell is sufficiently thin and the 
parameter u lies within the instability inter- 
val. The reasoning is as follows. 

Given the parameter u (= kyd/2), that 
meets the condition (j79|) . the thickness of the 
cell d, that enters Eqs. (f77 j) - (f7H|) through the 
parameters we = d/{2lg) and = d/{2l^), 
can be changed independently by varying the 
fluctuation wavenumber ky so as to keep the 
parameter u fixed. By this means the first 
two positive terms on the right hand sides of 
Eqs. (f77 |) -(f7S |) can be reduced to the limit 
where the sign of the total sum is determined 
by the last negative term. 

From the above discussion it follows that 
each value of u from the instability interval 
defines the critical point where the determi- 
nant of the critical mode vanishes. This point 
is characterized by the thickness of the cell 
and the fiuctuation wavelength, A = 2n /ky. 
The thickness can be computed as a func- 
tion of u by finding the positive root of the 
quadratic equation det Ma = 0. The wave- 
length then can be found from the relation 
u = kyd/2 = 7id/\ (see Eq. (gTj)). 

Geometrically, these calculations give a set 
of points {d, A) in the thickness-wavelength 
plane and by this means the instability inter- 
val turned out to be mapped onto the curve 

in the d-X plane. This curve represents 
the boundary of the instability region for the 
corresponding fiuctuation harmonics. 

It is now rather straightforward to carry 
out the described procedure and derive the 



parametric equations for the curve Fq, in the 
form 



D = Xa{u) = uXaiu)/^, 

A = A„(m), 



(80) 



+ [(1 - pYitanh^u + r^y 



I — ( 1 — p) (tanh^ u + 



Aa(M) = — -TT- 
+ [(1 - p)2(r^ tanh^M + 1)2 

where 



(81) 



( 1 — p) (r^ tanh^ u + 1) 



(82) 



ta{u) = q24.Pa{u){q24 " 1a{u)) , (83) 



If) h 



D is the size parameter; A is the wavelength 
parameter; and is the azimuthal anchoring 
parameter. 

In Eqs. dHHl)- (IHSl) the thickness and the 
wavelength are both conveniently scaled by 
the zenithal anchoring extrapolation length, 
Zg, and we shall use these formulas to per- 
form numerical analysis later on. But such 
analysis is not required for our subsequent 
discussion. 

Eqs. (jHn|) - (|HHjl describe the boundary of 
the instability region provided the instability 
interval is not empty. So, we turn back to 
the instability condition ()79|) to specify the 
regions of ^24 in which the instability may 
occur. To this end we discuss the behavior of 
the functions 7s (tt) and '~ia{u) that enter the 
right hand side of the inequality (ffn|) . The 
results, illustrated in Fig. El are not difficult 
to obtain analytically. 

Fig. |21^a) shows that, if the elastic con- 
stant K2 is smaller than Ki and r = 
K2/K1 < 1, the function 7^ increases from 
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7a(0) = 2r asymptotically approaching the 
value 7a (oo) = 4r/(r + 1) at large u, while 
7s is a decreasing function that varies from 
7,(0) = 2 to 7,(oo) = 4r/(r + 1). As 
is seen from Fig. |21^b), in the case where 
r > 1 the functions 7^ and 7, also monotoni- 
cally approach the common asymptotic value 
7a(cxD) = 7s(c)o) = 4r/(r + 1). But now we 
have 7,(0) = 2r > 7,(0) = 2. 

At the later stage we will see that there 
are two qualitatively different regimes of in- 
stability depending on whether the instability 
interval is bounded from above or not. So, we 
have the two characteristic values of ^24^ 



g(^) = min(2,2r), 



2(1 -p) 



4r 
r + 1' 



^5) 
^6) 



Now we pass on to discuss the stability of the 
planar structure in the intervals separated by 
these values. 



B. Stability diagrams 

In this section we study the stability of the 
planar structure by analyzing the behavior 
of the functions which define the boundary 
curve (jHUj) of the instability region. Qualita- 
tively, the analysis can be performed without 
resorting to numerical computations and we 
shall use the numerical results only for illus- 
trative purposes. 

Fig. 121 clearly indicates the interval 



< g24 < q, 



(1) 

c 



^7) 



as the region where the inequalities ()79p can- 
not be satisfied and the structure is stable. 
The stability conditions ()87|) are equivalent 
to the long known Ericksen inequalities Q 



derived in Ref . [2J| . Fig. |21 shows this stabil- 



ity region in the r-g24 plane. 
Next we consider the interval 

g«<g24<gf , (88) 

where, as is demonstrated in Fig. Et^a), the 
instability takes place at m < Wmax for the 



critical fluctuation mode which symmetry de- 
pends on the parameter r. For relatively 
small elastic constants K2 with r < 1, this 
mode is antisymmetric (see Fig. |21^a)). Re- 
ferring to Fig. Efb), it can be seen that in the 
opposite case with r > 1 the critical mode is 
symmetric. 

Since the instability interval is bounded, 
the size parameter D on the boundary 
curve (jHn|l is described by the function Xa{u) 
and reaches its maximum, Dc, as u varies 
from to Umax- It follows that in this regime 
the planar structure is unstable only in suffi- 
ciently thin films and dc {D^. = dc/le) is the 
critical thickness. On the curve this criti- 
cal point is also characterized by the critical 
fiuctuation wavelength Ac (Ac = Xc/h)- 

Since g24 = 7a at m = Umax, Eqs. dZZI) 
and (f7S|l show that the ^24 dependent term 
vanishes at u = and u = Wmax- When 
WgW^ ^ 0, the result is that Xa{0) = 
Xa{uma.x) = and the curve F^ forms the 
loop that encloses the region of instability. 
Fig. m^a) illustrates this result for r = 1.5 
when the boundary curve F, is determined 
by the symmetric critical fiuctuation mode. 

Now we consider the case of large g24 with 



924 > g?^- 



(89) 



In this case the instability interval is no 
longer bounded from above. From Fig. El it is 
clear that the instability condition (ffn|) does 
not impose any restrictions on u for the criti- 
cal mode, whereas the instability interval for 
the non-critical mode is: u > Mmin- Thus, the 
fiuctuation modes both lead to the instability 
at sufficiently large u. Eqs. flHT |) - (jH^ provide 
the limiting value of the functions Xs{u) and 
Xa{u) at M — > oo: 

Aoo = -{-(l-p)(r. + l) + [(l-p)' 
X (r^ + If + 4r^g24(g24 - )]'/'}, (90) 

so that the size parameter, D = mAq/tt, in- 
creases indefinitely in this limit. The result 
is that the planar structure is unstable at 
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any thickness of the cell and the instability is 
caused by the short wavelength fluctuations 
with A < Aoo- This result equally applies to 
the case of negative g24- 

It should be emphasized that using the 
continuum elastic theory we are limited to 
the case where the characteristic length of 
director distortions is large as compared to 
a microscopic scale of intermolecular dis- 
tances. This limitation can be taken into 
consideration by eliminating the fluctuation 
harmonics with large wavenumbers, k, 



x,y 



> 



/Cuv = 27r/Auv, from the Fourier series ex- 
pansion (jnZj), so that the wavenumber k^^ 
gives the large momentum (ultraviolet) cut- 
off in the theory. In this approach destabi- 
lizing effect of the critical short wavelength 
fluctuations may only be indicated if Aoo is 
larger than Auv Our results are scaled by the 
zenithal anchoring extrapolation length, Iq. 
Assuming that the anchoring is sufficiently 
weak and I0 ^ Auv, we conclude that in this 
case the effect of the cut-off on the accuracy 
of our results is negligible. 

Figs. Efb-c) illustrate transformations of 
the stability diagram in the D-A plane as 524 
increases beyond gc^^ at r > 1. The curve 
of the critical mode forms the boundary 
of the instability region similarly to the case 
of the interval ()88p. The curve of the non- 
critical mode Fa resides within this region 
and is not shown in the figures. Both curves 
rapidly approach the horizontal asymptote 
A = Aoo. 

We show in Fig. IH^b) that, when ^24 is 
close to gc , the shape of the curve F^ still 
somewhat resembles the loop developed in 
the interval (jHHj) . There are two points 
marked as -Dmin and -Dmax at which the 
curve bends and at which the function Xs{u) 
reaches its local minimum and maximum, 
correspondingly. Further increase of ^24 re- 
duces the distance between -Dmin and -Dmax 
and the points finally disappear after merg- 
ing at sufficiently large g24- Fig-IH^c) presents 
the stability diagram in this regime. 

In general, Fig.Elqualitatively summarizes 



the results of this section for the three differ- 
ent intervals of ^24 described in Eqs. (jH7|) - 
(jSni). The interval (|H^ was previously indi- 
cated as the instability region in Ref. 2l| and 
as the region where modulated structures can 
exist in Ref. @. 

Our results for the interval ()88|1 show that 
in this case the instability induced by the i^24 
term in cells of subcritical thickness may also 
lead to the formation of periodically mod- 
ulated structures. In addition, the critical 
wavelength Ac provides an estimate for the 
period of the emerging structure near the 
critical point. 



C. Effects of azimuthal anchoring 

The critical thickness and the critical 
wavelength are meaningful only when ^24 is 
within the interval (j88|l . By using Eqs. (j80|l - 
(l8Hj) they both can be computed numerically. 
The numerical procedure involves two steps: 

(a) solving the equation x'^{uc) = to find 
the maximum of Xa for the critical mode; and 

(b) evaluating = Xa{uc) and Ac = Aa(uc). 
We also present the results for Aoo computed 
from Eq. PH) . 

The critical thickness parameter, Dc, the 
critical wavelength parameter, Ac, and the 
parameter Aoo as functions of ^24 for differ- 
ent values of the azimuthal anchoring param- 
eter, = W^/Wg, are plotted in Figs. El El 
and m As it can be expected, the curves 
demonstrate that the structure becomes less 
stable as the azimuthal anchoring strength 
decreases. So, in order to estimate the criti- 
cal thickness from above, it is instructive to 
discuss the limit of weak azimuthal anchor- 
ing, W(f, 0. This limiting case was consid- 
ered in Refs. jl^ and we clarify some of 
the results by using our approach. 

It is, however, should be stressed that in 
this case the planar structure is marginally 
stable with respect to the long wavelength 
symmetric fluctuations with ky = regard- 
less of i^'24- Mathematically, it follows be- 
cause the determinant (|77|) vanishes at u = 
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and W(i, = 0. The reason is that all uniform 
planar structures are energetically equivalent 
in the absence of azimuthal anchoring. 

Nevertheless, we may apply our method 
to this case by taking the limit — > 
in Eqs. (|8(jp - (|9(Jj) . In this limit the expres- 
sions (|HT|). (jH^ and (jnOl) assume the simpli- 
fied form 



Xs{u) 



(1 — p) tanh^ m' 
27rtaiu) 



il-p) 
Aoo = AJoo) 



(91) 
(92) 
(93) 



From Eq. ()92|) it is seen that the behavior of 
the antisymmetric harmonics at m = does 
not differ from the case in which ^ and 
Aa(0) = Xa{0) = 0. But, for the symmetric 
mode, this is not the case. From Eq. (PT|) we 
have As — > cxD at u — > and 



= xs{0) = 2q. 



24iy24 
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(94) 



When r > 1 and the critical mode is sym- 
metric, Eq. provides the exact expres- 
sion for the critical thickness at = 0. In 
Refs. HQ the result ^ has been derived 
as an approximation. 

The stability diagrams in the D-A plane 
at r > 1 and = are shown in Fig. |H| In 
the opposite case of r < 1, so long as g24 < 2 
the diagrams are quite similar to those de- 
picted in Fig. m Otherwise, at g24 > 2, the 
non-critical symmetric mode will change the 
stability diagram presented in Fig. ^c). The 
result, however, is not too different from the 
diagram shown in Fig. IHl^c). 

If the azimuthal anchoring strength is not 
identically zero, Eq. ()94|) estimates the criti- 
cal thickness from above. In particular, sub- 
stituting into Eq. and taking the 
limit r ^ oo we arrive at the conclusion that 
the critical thickness cannot be larger than 
16le. 

Dependencies of the critical thickness, D^., 



and the critical wavenumber, 



2n/ Ac, on 



the ratio = W^/Wq for various values of r 
are plotted in Figs. IHl and ^| As is seen from 
Fig. ini the critical thickness declines steeply 
in the immediate vicinity of the origin. Typi- 
cally, the critical thickness at = appears 
to be halved at = 0.1, so that even for 
r = 20.0 we need very small azimuthal an- 
choring Vu, ^ 0.03 to have ~ 8 {dc ~ Sle). 
Referring to Fig. the critical wavenumber, 
Kc = 27r/Ac, also start growing rapidly, but 
the effect is less pronounced at large values 
of r. So, we have Kc ~ 0.2 {kc ~ 0-2 /Ig) at 
Tyj ^ 0.03 and r = 20.0. 

When the saddle-splay elastic constant 
meets a Cauchy relation: K24, = {Ki + 
K2)/2 [1, [in, we find the two intervals for 
the twist-splay ratio where the regime of in- 
stability defined in Eq. (jHHj) takes place: 3 < 
r < 3 + 2^/2 ^ 5.82 and 3 - 2^2 ^ 0.17 < 
r < 1/3. The critical thickness and the criti- 
cal wavelength as a function of the parameter 
r varying within these intervals are plotted in 
Fig.ini 

For r < 1/3, the critical fluctuation mode 
is antisymmetric and Fig. ITlT c) shows that 
the critical wavelength remains finite at = 
0. Referring to Figs. ITlT aVfb). the critical 
thickness at r < 1/3 is an order of magnitude 
smaller than in the case where r > 3 and the 
critical mode is symmetric. 

According to Ref. [2^, the latter presents 
the case in which the ratio r grows anoma- 
lously large due to an increase of the twist 
constant K2 in the vicinity of the nematic- 
smectic-y4 transition. In this case, the esti- 
mate of the critical thickness at = 0.0 

(2) 

and g24 = 5'c ~ 3.41 (r ^ 5.82) provides 
the upper bound for to be about 9.7/51. As 
is shown in Fig. ITlT b). this estimate can be 
significantly reduced in the presence of az- 
imuthal anchoring. 

D. Director fluctuations at substrates 
near the critical point 

In this section we apply the results of 
Sec. IIIIBI to study the correlation functions of 
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director fluctuations at the plates bounding 
the cell. Speciflcally, we shall use Eqs. (f7^ - 
()75|) that express the correlator in terms the 
inverse of the matrices and given in 
Eqs. dH) and (I59E1) . 

Obviously, the results obtained in the 
Gaussian approximation are inapplicable in 
the region where the structure is unstable. 
So, we shall restrict our considerations to the 
stability region and study what happen when 
either the thickness or g24 varies so as to get 
closer to the boundary of the instability re- 
gion. 

We start with the discussion of the long 
wavelength limit, u 0, for the matrices 
M~^. The result of calculations is that these 
matrices are diagonal at m = 0: Ad~^ = 
diag(w^\ + and M~'^ = diag{{w^ + 
r)^^,Wg^). It is seen that in this limit the 
correlator diverges for vanishing the anchor- 
ing strengths, W^f, = or Wg = 0. This is a 
consequence of the marginal instability dis- 
cussed in the preceding section. 

As far as the large wavenumber (short 
wavelength) limit is concerned, it can be 
shown that in the stability region the matri- 
ces (and the correlator) both decay to zero 
at M oo. Interestingly, this is not the 
case at the boundary of the instability in- 

(2) 

terval when ^24 = or g24 = Qc ■ In this 
case we have the non-zero limits M~^{oo) = 

{wg + w^y^ ^2]^ ^"""^ for g24 = and g24 = 

qi^\ respectively. This anomaly is a precur- 
sor of the instability caused by short wave- 
length fluctuations. As a result, using the 
widespread approximation with K24 = does 
not give the correlation functions that prop- 
erly behave at large wavenumbers. 

We demonstrate in Figs. IT^ and IT^ which 
show the spectra of the critical and the non- 
critical fluctuation modes computed from 
Eq. (|75|) at r = 1.5, that the critical increase 
of the symmetric fluctuation mode becomes 
sharply peaked at the critical wavenumber as 
the thickness of the cell approaches its criti- 
cal value. In addition. Fig. IT^ shows that the 



symmetric and the antisymmetric modes are 
dominated by the in-plane and by the out-of- 
plane fluctuations, respectively. 

V. DISCUSSION AND CONCLU- 
SIONS 

In this paper stability of the uniform di- 
rector distribution in a planar NLC cell has 
been studied in the presence of the saddle- 
splay term. We have shown that the phys- 
ical mechanism behind instabilities induced 
by the ^^^24 term is governed solely by the di- 
rector fluctuations at the substrates that can 
generally be characterized by the free energy 
and the probability distribution. We have de- 
veloped the approach to study such instabil- 
ities by using the correlation function of the 
fluctuation fleld induced by the fluctuations 
at the conflning wall (the surface part of the 
correlator) for derivation of the stability con- 
ditions. 

This approach in combination with the 
mirror symmetry considerations has been ap- 
plied to the case of NLC planar cell. It is 
found that there are two types of fluctuation 
modes, that we have called symmetric and 
antisymmetric modes, depending on the par- 
ity under the mirror symmetry transforma- 
tion involving reflection in the middle plane 
of the cell. 

We have devised the analytical method 
to analyze the stability conditions for the 
fluctuation modes of different symmetry. In 
this method the thickness of the cell and the 
fluctuation wavelength form the plane and 
the boundary of the instability region in this 
plane is described as a curve deflned in the 
parametric form. It appears that no numer- 
ical calculations is required as far as quali- 
tative analysis of stability diagrams is con- 
cerned. The analysis revealed the two dif- 
ferent regimes for instabilities caused by the 
saddle-splay term depending on the value of 
K24. 

When ^^^24 falls within the range between 
2mm{Ki,K2) and 4KiK2/{Ki + K2), the 
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planar orientational structure loses its sta- 
bility only in sufficiently thin cells and the 
critical thickness, dc-, together with the wave- 
length of the critical fluctuation mode, Ac, 
determine the critical point. In this case, as 
is shown in Fig. [121 the spectrum of critical 
fluctuations at the substrates grows sharply 
peaked at the critical wavelength when ap- 
proaching the critical point. The period of 
modulated structure emerging at the critical 
point is thus determined by the critical wave- 
length. 

This K2A induced instability takes place at 
any elastic anisotropy parameter r = K2/ Ki. 
The sole exception is the case of elastic 
isotropy in which Ki = K2. In contrast with 
the periodic s play -twist Freedericsz transi- 
tion w 
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where spatially modu- 
lated pattern comes into play only at suffi- 
ciently small twist-splay ratio with K2/K1, 
the surface elasticity driven instability can- 
not be hindered by large elastic anisotropy, 
but rather, as is shown in Fig. IHl the greater 
elastic anisotropy the larger the critical thick- 
ness can be. The symmetry of the critical 
fluctuation mode, however, depends on the 
parameter r: the mode is antisymmetric at 
r < 1 and is symmetric in the opposite case 
of r > 1. Fig. ^1 shows that the in-plane 
and out-of-plane fluctuations prevail depend- 
ing on the symmetry of the fluctuation mode. 

The azimuthal anchoring turned out to 
have a profound effect on both the critical 
thickness and the critical wavelength. The 
absence of the azimuthal anchoring presents 
the limiting case where the planar structure is 
marginally unstable with respect to the long 
wavelength fluctuations with ky = regard- 
less of the K24: term. As a consequence, the 
critical wavelength Ac increases indefinitely 
in the limit of weak azimuthal anchoring, 
Ac — 00 at — > 0, while the limiting value 
of the critical thickness can be computed ex- 
actly (see Eq. (jHH)) giving the upper bound 
for the critical thickness. 

In Sec. lIVCJl the absolute upper bound for 
the critical thickness was found to be 16 Iq. 



Fig. ini shows that the critical thickness could 
have been significantly reduced in the pres- 
ence of relatively small amount of the az- 
imuthal anchoring. 

Using the Cauchy relation: K24, = {Ki + 
-^2)72 ja, we have found that the in- 
stability may occur at both sufficiently low 
and high twist-splay ratios with r < 1/3 and 
r > 3, respectively (see Fig. [TT] and the re- 
lated discussion at the end of Sec. II V C|) . 

The corresponding condition for the peri- 
odic Freedericsz transition requires the ratio 
r to be below Tc ~ 0.303 [4fl,|41i] which places 
a slightly more stringent constraint on the 
value of r than the above inequality: r < 1/3. 
But, as is seen from Fig. ITlT a). for small r, 
the critical thickness is an order of magni- 
tude smaller than in the case of large r pre- 
sented in Fig. llir b). where dc can be of order 
of several microns provided the extrapolation 
length I0 varies in the range 0.1 — 1 fim. 

For typical nematics, the twist-splay ra- 
tio does not exceed unity, but close to the 
nematic-smectic-y4 transition the parameter 
r becomes anomalously large y| leading to 
the K24 induced instability of the ground 
state at r > 3. This may result in the appear- 
ance of modulated orientational structures as 
suggested in Ref. j23l |. 

When 7^24 > 4KiK2/{Ki + K2) or i^24 < 
0, the short wavelength fluctuations with A < 
Aoo will render the planar structure unstable 
at any thickness of the cell. With the Cauchy 
relation such instability will take place when 
the twist-splay ratio is either less than 0.17 or 
greater than 5.82. In contrast with the above 
discussed regime, this instability though does 
not impose any restrictions on the ffim thick- 
ness, in general, cannot be unambiguously re- 
lated to the periodic pattern formation. This 
case requires a more detailed additional study 
of orientational structures in the instability 
region where the Gaussian approximation is 
inapplicable. 

Our concluding remark concerns the gen- 
eral method for separating out the contri- 
bution of director fluctuations at confining 
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walls to static correlation functions. We have 
demonstrated that this method can be used 
as a useful tool for studying orientational in- 
stabilities in confined liquid crystals. For 
this purpose, we have restricted ourselves 
to the case of uniaxial director fluctuations 
with uniformly distributed degree of order- 
ing. But a complete treatment of fluctuations 
at confining surfaces is required in studies 
of such phenomena as electrohydrodynami- 
cal pattern formation i^l, instabilities un- 
der shear flow (4, , wetting and back- 
flow 1^. These more general considerations 
involving fluctuations of spatially varying or- 
der parameter tensor coupled to the transla- 
tional degrees of freedom is well beyond the 
scope of this paper and we will extend on this 
subject elsewhere. 



ergy of the fluctuation harmonics in the fol- 
lowing matrix form 



d/2 



d/2 



^l=±l 



z=fid/2 



(A3) 



Using the transformation (j39|l will render all 
matrices in Eq. ()A3j) real-valued. These are 
given by 



Bi = -ky 



K2 
Ki^ 

Ki 
7 ' 



(A4) 
(A5) 
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APPENDIX A: FREE ENERGY OF 
FLUCTUATIONS IN NEMATIC CELL 

In this appendix we briefly comment on 
the derivation of the expression ()4()j) for the 
free energy. We start with the second order 
expression for the free energy (jH)): 



1 



" 2 jy 



(Al) 



i?(2) 



1 



2 Js^ 

K2A [Ody 



^dyO] 



ds, 



(A2) 



where we retain the director derivatives with 
respect to x for discussing its role later on. 

Substituting the expansion (jHTf) into 
Eqs. f)Al|) and ()A2|) will provide the free en- 



Kikl 







D. 



+ 



K2kl 



fXky 



K.kli. 



(A6) 
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w. 



(m) 







w, 



(m) 



(A7) 



where wf and W^g'"'^ are the azimuthal and 
the zenithal anchoring strengths at the sub- 
strate z = fid/2. 

For a symmetric cell with identical sub- 



\8trates, we have W. 



(±) 



and W, 



(±) 



Wo 



This is the case we deal with in this paper. 
But the generalized expression ()A7|) can be 
used to discuss the effects of anchoring en- 
ergy asymmetry. 

Another generalization is that the deriva- 
tives with respect to x have not been ne- 
glected from the very beginning and we have 
the term proportional to /c^ in the expression 
for the matrix ()A6|1 . 

But the stabilizing magnetic field H = 
He^ likewise will change the matrix ()A6|1 : 
Ci — > Ci+XaH"^!, where Xa is the anisotropic 
part of the magnetic susceptibility. It follows 
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that the effects of non-vanishing wavenum- where 
bers kx and of the stabihzing magnetic field 
are equivalent, so that, at 7^ 0, the struc- 
ture is more stable than in the case with 
kj. = 0. Thus, the fluctuation harmonics with 
kx ^ do not affect the results of stability 
analysis and can be safely eliminated from 
the consideration. 

Assuming that k^ = and the cell is sym- 
metric, Eqs. flA3|) - ()A7|) can now be rewritten 
in terms of the dimensionless parameters de- 
fined in Eqs. and (jTTj) to yield the ex- 
pressions fpH|) -(|i ^ . 



APPENDIX B: GREEN FUNCTION 
METHOD FOR NEMATIC CELL 

In this appendix we deduce an expression 
for the correlator of director fluctuations in 
the NLC cell by using the Green function 
formalism described in Sec. Ill Bl We shall 
find that this method can also be employed 
to rederive the result for the surface part of 
the correlator given in Eqs. (|68|) and (|69p. In 
addition, this appendix provides a number of 
technical details omitted in the bulk of the 
paper. 

We start with the expression for the free 
energy functional taken, similarly to 

Eq. ()10|) . in the following form 



-u 



ip^Ltl^ dr 



fi=±i 



r=fiu 



(Bl) 



(B2) 



Substituting the fluctuation field (jH^ ex- 
pressed in terms of the matrices that 
satisfy the Euler-Lagrange equations (jiH|) 
and meet the boundary conditions (j63p into 
Eqs. ()B1|) and (j42j) gives the matrix that 
enters Eq. (jUBI) 



A^ 



N+- 



(B3) 



A^, 



fii/ 



(B4) 



(B5) 



=flU 



We define the Green function as the solu- 
tion of the following boundary value problem: 



LGiT,T') = -6ir-r')i, 



Q±G 



0. 



--±u 



(B6a) 
(B6b) 



Since the operator K from Eq. ()23a|) is 
proportional to the operator —L: K = 
—uKiS/d L, we need to modify the rela- 
tion linking the Green function and the 
correlator as follows 

Next we seek the Green function G in the 
form of a sum 

G{T,r') = G'^'\T,r')+G^^\T,r'), (B8) 

where is the Green function of the 

Dirichlet problem: 



LG'^'\t, r') 
G^''\±u,t'\ 



r')L 



0. 



(B9a) 
(B9b) 



The boundary conditions ()B9b|) of the Dirich- 
let problem present the limiting case of strong 
anchoring and the Green function is pro- 
portional to the bulk part of the correla- 
tor (lB7|l . 

Since the sum ()B8|) is the solution of the 
problem ()B6|) . the equations and the bound- 
ary conditions for G^^^ are given by 



LGW = 0, 



T=fJ,U 



-uQfG^'^ 



(BlOa) 
(BlOb) 



T=fJ,U 
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Clearly, the surface part of the Green 
function represented by G^^^ accounts for 
the difference between the boundary condi- 
tions ()B6b|) and the strong anchoring condi- 
tions (lB9E|l . 

We have already pointed out the difference 
between the free energy (|iup and the func- 
tional ()B1|) above the modified relation ()B7|) . 
Likewise, this difference will slightly change 
the Green formula (fT^ 



\il)~^Lip — (p^Lt/) dr 

E k^^^V - _ • (Bii) 



and compute the matrices Q^'^\t') from the 
boundary conditions ()B10b|) . The right hand 
side of Eq. flB10b|) can be simplified by means 
of the relation ()B12|) . Then substituting 
Eq. (pT5|l into Eq. ()B10bl) and using Eq. (|B4|l 
yield the equations for Q^'^\t') in the final 
form 



N,.^^''\t') = u[¥^\r't . (B16) 

u=±l 



At=±l 



It can now 
Eqs. dHU), 



be 



easily verified 
and (fRT6|l 



that 



There are two important relations that al 
most immediately follow from the Green for- 
mula dEIH): 



T, T 



()B15;i ana ifmhji pro- 
vide the expression for the surface part of 
the correlator given in Eqs. ()68|) and ()69p. 

The relation (IbT21) can also be applied to 
derive the bulk part of the Green function 
by means of the so-called Wronski construc- 
[¥^'\T')f (B12) tion 0. The result is 



-flU 



and 



(B13) 



Eq. (jB12|l can be derived from Eq. (jBllj) by 
setting V = G^^\t,t') and y = ^(^)(r)c^('^), 
whereas for the relation ()B13|) the substi- 
tution is: tl) = ^(''^(r)?/?^'') and (fi = 



T)cp 



(p) 



Combining Eqs. ()B3|) - ()B5|) and the rela- 
tion ()B13jl shows that the matrix is sym- 
metric 



R = 



^(-)(r)i?+[^(+)(r')]" 



-[N. 



T<r' 
t>t'' 
(B17) 



Our final remark concerns the symmetry 
relation 



P¥''H-t)P 



(B18) 



N. 



(B14) 



So, the matrices that enter the right hand 
sides of Eqs. (IB^ and (jU^ are identical. In 
addition, the matrices Mq, defined in Eq. (jSHl) 
are also symmetric. The latter, however, can 
be seen from the explicit formulas ()59|) . 

Now we write the surface part of the Green 
function as the general solution of Eq. (jBlOajl 

G(«)(r,r') = J2 ^^"\r)e^"\r') (B15) 
i/=±i 



that follows from Eqs. 
combination with Eqs 
the following result 



N 



(b) 



aj) and ()64b|) and in 
Bill and (jHH) gives 



(B19a) 
(B19b) 



This result simplifies algebraic verification of 
the relation (|U7|). 
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